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ABSTRACT 

The investigation dealt with the problem of biorhythmic 
criticality and its influence on human error and accidents, 
based on data from 4346 naval aircraft mishaps in the Fiscal 
Years 1968-1973. Observed occurrences of mishaps were com- 
puted under different aspects and compared against expected 
occurrences, obtained from a mathematical model. 

The findings of the study supported the relativity- idea 
in the theory of Biorhythms: critical states of the three 

biorhythmic cycles have different effects on groups of people 
with different characteristics. For pilots, the study re- 
vealed no significant influence from Biorhythms by straight- 
forward application of critical days, as observed in previous 
research. Significant results at varying levels between a = 
0.15 and 0.03 were obtained when considering different age 
groups among pilots. The intellectual cycle was found to be 
irrelevant in connection with aircraft mishaps. A high nega- 
tive correlation was discovered concerning the deviation be- 
tween observed and expected mishaps, when pilots under thirty 
years of age were compared with those above thirty (a = 0.117). 
The question of increased significance when applying critical 
intervals rather than critical days as defined in the concept 
of the "Critical Category" was answered negatively. 

Final goal for the application of Biorhythms in that field 
is considered to obtain a predictive device for operational 
purposes, reducing the accident rate by avoiding flying days 

for pilots with disadvantageous criticality states. 
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I. INTRODUCTION 



A. THE PROBLEM 

For a long time, since the beginning of industrialization 
in the Nineteenth Century, it was well accepted that the bot- 
tleneck of any system was the machine: high costs, limited 

availability, technological limitations, and unsatisfactory 
reliability of its parts and therefore the whole machine con- 
stituted severe constraints in the productivity of any system. 
It was only natural that all concentrated efforts of the in- 
dustrialized world were devoted to improvements of this "weak 
link". Mankind proved to be successful in that respect; the 
only uncomfortable experience, however, was that the man him- 
self had taken over the role of the weak link as machines be- 
came better and more complex. People also began to realize 
that it was a much more complex and sophisticated problem to 
remove the bottleneck "man" than it had been in the machine's 
case. What influences his performance in the necessary man- 
machine interaction? How can one predict his behavior and 
performance variation in trying to accomplish his task? These 
were some of the questions man presented in the man -machine 
interface . 

Another drastic change occurred in our attitude toward 
'man' and 'machine'. It was only about one hundred years ago 
that loss of a machine in many parts of the world was con- 
sidered to be much more severe than loss of a man (worker). 
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Today we are not only concerned about good performance of the 
man but also about his safety and his welfare as a human 
being. Considerations like these led to the question: Being 

unable to eliminate it, is there any possibility to predict 
at least the variations in human performance? If we could 
predict times of 'high' or 'low' performance, relative to a 
certain (possibly unknown) mean-level, many aspects of human 
life would be positively affected, making use of this know- 
ledge. One of the numerous attempts to solve this problem 
was the concept of Biorhythms, first brought up at the end 
of the Nineteenth Century in Germany and Austria. 

B. THE THEORY OF BIORHYTHMS 

The concept of Biorhythms claims that each human being is 
influenced throughout his life by three different cycles, 
having a different period: 

A 23-day period cycle, governing the physical 

condition . 

A 28-day period cycle, governing the emotional 

condition. 

A 33-day period cycle, governing the intellectual 

condition. 

All three cycles can be imagined to have the shape of a sine- 
wave, having a positive and a negative half -cycle. They all 
start simultaneously at birth, moving upwards, and continue 
unchanged through man's life. According to research results 
in various fields, accidents occur more often on those days 
where a cycle starts over again, or crosses the "zero- line" 
on the way from high to low |Ault and Kinkade, 1972; Thommcn , 
1973; Senzaburo, 1969]. These days are called "critical days" 



9 



Consequently, there are three kinds of critical days, depending 
on how many of the cycles intersect on the axis on a given day. 
In addition to the already mentioned critical day (one cycle 
intersects), there is a "double-critical day" (two cycles in- 
tersect) and a "triple-critical day" (all three cycles inter- 
sect) . Even higher accident possibility is predicted for the 
multiple critical days, especially those involving physical- 
emotional criticality [Willis, 1972], 

physical 

emotional 

intellectual 

© © critical days 




Figure 1. The Three Biorhythmic Cycles (Starting at Birth). 

The 23-day physical rhythm has a positive and a negative half 
cycle of 11.5 days each, indicating better physical condition 
in the plus period and reduced physical capability in the 
minus -period . The emotional cycle with a period of 28 days 
and two 14 day long half cycles, predict optimism and posi- 
tively influenced emotions during the plus period, emotional 
instability and pessimistic attitudes during the minus -period. 
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The intellectual cycle with a period of 33 days and two half 
cycles of 16.5 days length each are considered to be of minor 
importance to accidents, having more influence on mental 
alertness in problem solving and similar tasks. 

It must be emphasized that biorhythmic theory does not 
try to predict absolute "states", it gives only relative 
classifications like a nominal scale does. That means that 
one person on a critical day, which is considered to be very 
disadvantageous still could perform much better in a certain 
task than another person, having the respective cycle at its 
maximum, just because of different initial capability. This 
could very well lead to a point that persons with above aver- 
age capabilities in reaction time, motor response, etc. are 
able to avoid an accident in a critical situation, where the 
majority would have suffered an accident. This point will 
become important later in the study. 

C. HISTORIC BACKGROUND 

History books tell us that about 400 B.C. Hippocrates 
advised his students to observe "good" and "bad" days among 
the healthy and ill, and to take those fluctuations into con- 
sideration in the treatment of the patients. 

Between 1897 and 1902, psychologist Dr. Hermann Swoboda 
(University of Vienna) did some initial research on recurrence 
of pain and swelling of tissues. Medical review led to the 
discovery of the 23-day and the 28-day cycle. Also, the first 
biorhythmic slide rule for determination of the critical days 
was designed by him. 
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The nose and throat specialist, Dr. Wilhelm Fliess, M.D., 
discovered independent from Swoboda these two cycles, based 
on diagnosis and observations of patients. He called the 
23-day cycle the 'masculine', the 28-day the 'feminine' cycle. 

The third cycle (intellectual, 33 days) was discovered 
much later in the 1920 's by the engineer, Dr. Alfred Teltscher 
(University of Innsbruck, Austria). He based his findings on 
academic performance fluctuations of high school and college 
students, caused by periodic secretions of glands affecting 
the brain cells. The same result was obtained by Dr. Rexford 
Hersey (University of Pennsylvania) in a study using data 
from workers in railroad shops. 

About the same time people began to think about biorhythmic 
applications in the field of accident research. The mathema- 
ticians Dr. Alfred Judt (Bremen, Germany) and H. R. Frueh 
(Switzerland) had provided the necessary computational tools 
in form of calculation tables and hand-operated calculators 
for determination of biorhythmic criticality. In 1939, the 
first intensive study with a data base of 700 accidents from 
insurance companies was performed by Hans Schwing, a student 
at the Swiss Federal Institute of Technology in Zurich, 
Switzerland [Schwing, 1939], Schwing came up with the follow- 
ing results: Taking into consideration only the physical and 

the emotional cycle, he found 401 accidents occurring on crit- 
ical days, which is about 57%. Of these, 322 fell on single 
critical days, 74 on double critical days, and 5 on triple 
critical days. A similar strong significance of biorhythmic 
criticality on accidents was reported by Rcinhold Bochow and 
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Dr. J. Sennewald (Humboldt University, Berlin, [Bochow, 1954]) 
investigating 497 accidents of workers using agricultural 
machinery. Of these, 26.61 of the accidents occurred on 
single critical days, 46.5% on double critical days, 24.7% 
on triple critical days, and only 2.2% on non-critical days. 
The only known accident investigation in the field of avia- 
tion based on biorhythmic theory was performed with a data 

base obtained from the Guggenheim Aviation Safety Center, 

) 

Cornell University, involving private pilots. Of the ob- 
served accidents, 80% occurred on critical days of the pilot 
[Thommen, 1973]. 

In the past twenty years, numerous private and government 
institutions and businesses have made use of the Theory of 
Biorhythms in connection with accident prevention, especially 
in Japan, Switzerland, and West Germany. Nothing is known 
about research and application of this theory behind the iron 
curtain, except that extensive literature fell into Russian 
hands after World War II, including many of the findings of 
Dr. Swoboda [Thommen, 1973]. In recent years, several studies 
at universities in the United States have been reported 
[Willis, 1972]. 

Based on these reported results it was felt that the pos- 
sibility of a significant reduction of accident rates in var- 
ious military fields justified a thorough research in the 
matter . 
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II. PURPOSE 



Accident prevention is an important aspect everywhere in 
the military. Life, welfare, and health of the people are 
involved as well as the fact that accidents cause high (and 
possibly unnecessary) costs to the budget. This research 
investigates the possibility of using biorhythmic criticality 
for prevention of aircraft mishaps, based on 4,346 aircraft- 
mishaps which occurred in FY 1968-1973 in the U.S. Navy. It 
was felt that strong emphasis should be given to the follow- 
ing statement: The research investigated the usefulness and 

applicability of the biorhythmic concept, as it was found to 
be applicable by the many scientists and institutions or 
businesses, mentioned in the Historical Background. It did 
not try to find reasons and causalities which lead to the 
phenomenon of biorhythms. Therefore, only the value of it as 
a predictive tool for people's accident- likelihood as a func- 
tion of the day was under question. 

One of the biggest problems previously in investigating 
the applicability of biorhythms was the lack of appropriate 
data; also, it was just about twenty years ago that the 
availability of the computer and appropriate statistical tech- 
niques made it possible to look at sufficiently large sample 
sizes to answer the question: Do accidents occur by chance, 

or do they occur (significantly different from chance) ac- 
cording to biorhythmic criticality? The purpose of this 
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research was to answer this question for the specifi 
of aircraft-mishaps and the involved pilots. 



category 
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III. THE MATHEMATICAL MODEL 



Assume an aircraft mishap has occurred on a given day. 
Then, the criticality of this day in terms of the biorhythmic 
condition of the involved pilot can take on eight different 
mutually exclusive and qualitatively defined values: 



NC = non-critical 
P = physical critical 
E = emotional critical 
I = intellectual critical 
PE = double critical physical/emotional’ 

PI = double critical physical/intellectual 
El = double critical emotional/intellectual 
PEI = triple critical physical/emotional/intellectual. 

According to that, define a random vector X. = (x, . x_ 

-J K lj 2j 

x 8j) = (X i_j) describing the biorhythmic criticality of a 

given mishap day j , where 



1 




NC 

P 

E 

I 

PE 

PI 

El 

PEI 



and 



16 



1 if the mishap-day j has criticality i 



x 




0 



if the mishap-day 



has criticality other than i. 



\ 

Then it is possible to compute the probabilities p^, which 
are associated with the occurrence of a mishap-day j having 
criticality i, under the null-hypothesis that aircraft -mi shaps 
occurred by chance and are not influenced by biorhythmic crit- 
icality, such that p^ = P(x ) . 

Assume that in general, a set of n mishaps in form of 
appropriate data is available (i.e., giving birthday of the 
pilot and day of the mishap). Then, let 



n 

x. = E x. . , 

1 j = l ^ 

and 

n 

X = E^ X. = (x. ) = (x-j^ x 2 x g ) , j = 1, ,n. 

The random vector X is distributed according to the multi- 
nomial distribution with parameters n and (p^ p ? Pg) , 

and the following properties: 



RANGE Ry = {(x^....Xg): x^ = 0,1,. ...,n; i = 1,2,.. ..,8; 

8 

E x. = n } . 
i-1 1 



PROBABILITY MASS FUNCTION 



P x (x) 




x 



1 



X 



2 



*P 



1 



P 2 



x 



8 



^8 
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n! 



8 

n x. : 
i=l 1 



1 2 
Pi *P2 ' 



for 



x e B x , 



0 5 Pi 5 1> 

8 

Z 

i=l 



2 p. = 1, 



8 

Z 

i= 1 



Z x. - n, 

i 



MOMENTS: E [X ± ] = n*p i ; V [X ± ] = np^ 

and E [X] = n(p 1 p 2 p g ) . 

This model will be referred to as the "base model" later in 
the study. The distinction seemed to be necessary because 
several slight modifications in the dimensionality of the 
vector X had to be applied in order to investigate problems 
involving differently defined categories i 1 . The modifica- 
tion consisted in these cases of combining appropriate cate- 
gories i together, thus reducing to less than eight dimensions, 
but maintaining the distributional characteristics of the 
model. In general, it always is based on a multinomial dis- 
tribution of the kind 



n 



P x (x) - 



X 1 X k 



P- 



• P- 



A k 

Pk ’ 
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whe re 



x e i?x , 

o 5 Pj 5 i; 



k 

z 

i= 1 



1 . 



Whenever this modification has been used, a short explanation 
about it and the necessary new definitions of the involved 
parameters have been added. As an example, the random vector 
Y = (Y^ Y 2 Yj Y^) describes the criticality of a mishap day 
in terms of the four mutually exclusive categories: non- 

critical, single critical, double critical, and triple 
critical . 

Accordingly, p i = PCY^), i = 1, ,4; j = 1, ,n; 

represents the probability that a mishap day j has criticality 
i, following the same rationale as the base case, except that 
now Y = f(X) 3 Y 1 = X 1 

4 

Y, = Z X. 

2 i-i 1 

7 

Y _ = Z X. 

3 • r x 
1 = 5 

Y 4 ' X 8' 



A. PROBABILISTIC ASPECTS 



As already mentioned, the distribution of the random vec- 
tor X (under the Null -Hypothesis of mishaps occurring by chance 
and unaffected by biorhythms) was multinomial with Probability 
vector p = (p. .... p 0 ) . Therefore, the probabilities p- = 

i. O -» 
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P(X^j) had to be computed to use this null-distribution in 
the further statistical analysis. A set theoretic model 
was employed to obtain this information. Let 



A 


= set 


of 


all 


B 


= set 


of 


all 


C 


= set 


of 


all 


D 


= set 


of 


all 


Then 








BOC 


= set 


of 


all 


bdd 


= set 


of 


all 


COD 


= set 


of 


all 


BHCflD 


= set 


of 


all 


The ba 


sic appli 


cat 



non-critical days; 
physical critical days; 
emotional critical days; 
intellectual critical days. 



double critical 
double critical 
double critical 
triple critical 
ion of Venn Diag 



days, physical/emotional; 
days, physical/intellectual; 
days, emotional/intellectual 
days . 

rams might illustrate this. 





Figure 2. 



Set-Theoretic Model for the Calculation of the 
Probability-Vector for the Null-Distribution. 
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Employing the necessary set theoretic relationships, the 



elements of the probability vector £ were obtained as follows 



p^ = P (non-critical) 

= 1 - [P (B) + P(C) + P (D) - P(BOC) - P(BflD) - P(COD) 

+ P(BOCnD)]. 

P 2 = P(single critical physical) 

= P(B) - P(BOC) - P (BAD) + P(BncnD). 

p^ = P(single critical emotional) 

= P(C) - P(BHC) - P(CHD) + P(BOCnD). 

p^ = P(single critical intellectual) 

= P (D) - P(BflD) - P (COD) + P(BnC0D). 

Pi- = P(double critical physical/emotional) 

= P(BOC) - P(BOCnD). 

p^- = P (double critical physical/intellectual) 

= P(BHD) - P(BOCAD). 

p^ = P(double critical emotional/intellectual) 

= P (OOD) - P(BnCDD). 

Pg = P(triple critical day) = P(B/>CAD). 

The computation of the numerical values of the p_^ was based 
on the following rationale (explained for an example of the 
23-day cycle, the others are analogous). 

The set B corresponds to all critical days of the 23-day 
physical cycle. The number of days until the occurence (ex- 
clusively) of the next critical day of this type, given it 
has just occurred, is exactly the length of one half-cycle, 
namely 11.5 days. Therefore, one out of the 11.5 clays in a 
physical half-cycle is a critical one. Under the assumption 



21 



of the Null-Hypothesis that there is no influence from bio- 
rhythmic criticality, an accident (having occurred) is as 
likely to have occurred at a critical day as well as on any 
one of the other days in the half-cycle, thus leading to the 
idea of a uniform distribution: the probability that the 

accident occurred on the critical day is just 1/11.5 because 
the life of any man can be thought of consisting of a large 
number of these half- cycles , the numerical value will remain 
unchanged, thus yielding the wanted probability that an ac- 
cident occurs on a physical critical day by chance. 



Event 


Set 


Length of 
Half -Cycle (1) 


Probabi 1 ity (Set) = 1/ (1) 


Physical 


B 


11.5 


P (B) - * = 0.08696 


Emotional 


C 


14 


P(C) = = 0.07143 


Intellect. 


D 


16.5 


P (D) = jg-g- = 0.06061 


Phys/Emot 


BrtC 


11.5 x 14 = 
161 


P(BflC) = 1 ^ 1 = 0.00621 


Phys/Int 


BAD 


11.5 x 16.5 = 
189.75 


p ( B A D ) 189.75 

0.00527 


Emot/Int 


CAD 


14 x 16.5 = 231 


P(CflD) - * - 0.00433 


Phys/Emot/ 

Int 


BfiCfyD 


11.5 x 14 x 

16.5 = 2656.5 


P(HACnD) = 26 l 65 = 
0.00038 



Table I. Computation of Probabilities for the Nul 1 -Distribution . 
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Using these results the numerical values of the unknown prob- 
ability vector could now be computed: 



As it should be, = 1. Also, the probability that a mis 

hap falls on a critical day of any type whatsoever is 1 - p^ 
0.20356, which is in accordance with the value used in the 
cited literature and previous research [Thommen, 1973] 

B. COLLECTION AND PROCESSING OF THE DATA 

The data base consisted of 4346 aircraft mishaps, having 
been recorded in the fiscal years 1968 through 1973 in the 
files of the Naval Safety Center, Norfolk, Virginia. The dat 
contained pilot's name, cause code, date of birth, date of 
mishap, and was stored on magnetic tape. Because of the size 
of the data set, it would have been impractical and too time 
consuming to compute the criticality of a mishap day by hand. 
Therefore, a computer program had to be written to provide 
this information. The basic idea was to start with the birth 
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date of the pilot, which caused the j ' th mishap, and to find 
the state of the three cycles at the mishap day. This yielded 
three numbers, which were compared with the pre-defined crit- 
icality definitions, using logical I F- statements . According 
to this procedure, one of the possible eight criticalities 
was found to fit and printed behind the name of the pilot 
(such as "P" for a single critical day physical; if the mis- 
hap occurred on a non-critical day, the space provided for 
this information was left blank) . A simple counter then sum- 
marized the number of occurrences of each criticality category 
for the 4346 mishaps. 

Presentation of the computer program itself here was not 
considered to be of importance because of its triviality and 
the fact that a large portion of it just consisted of the 
necessary programming procedures for extraction of the rele- 
vant data portion from the master tape. More illustrative 
should be a sample of the actual printout, showing the display 
of the desired information. 
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Figure 3. Computer Output Sample 



